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ABSTRACT 


During  the  AMPTE  (Active  Magnetospheric  particle  Tracer 
jixplorers)  solar  wind  lithium  release  on  September  11,  and  again  on 
September  20,  1984,  an  intense  burst  of  electrostatic  noise  was 
observed  near  the  upstream  edge  of  the  ion  cloud.  Comparisons  with 
measurements  by  the  IMP-6  and  ISEE-1  spacecraft  show  that  the  spec¬ 
trum  and  overall  features  of  this  noise  are  very  similar  to  electro¬ 
static  noise  observed  at  the  earth’s  bow  shock.  A  stability  analysis 
using  realistic  parameters  shows  that  the  electrostatic  noise  can  be 
accounted  for  by  an  ion  beam-plasma  instability  caused  by  the  solar 
wind  proton  beam  streaming  through  the  nearly  stationary  lithium 
cloud.  The  growth  rate  of  this  instability  is  largest  when  the  ion 
density  and  solar  wind  proton  density  are  similar,  which  explains  why 
the  noise  only  occurs  near  the  outer  edge  of  the  ion  cloud. 

The  similarity  to  the  noise  in  the  earth's  bow  shock  suggests 
that  a  shock  me  exist  in  the  solar  wind  plasma  flow  upstream  of  the 
ion  cloud.  If  the  noise  is  associated  with  a  shock,  then  it  must  be 
an  electrostatic  shock,  since  the  ion  cyclotron  radii  are  too  small 
for  the  existence  of  a  MHD  shock.  Since  the  electrostatic  insta¬ 
bility  occurs  at  phase  velocities  near  the  lithium  thermal  velocity, 
the  electrostatic  turbulence  may  play  a  role  in  heating  the  lithium 
ions  and  transferring  momentum  from  the  solar  wind  to  the  ion  cloud. 
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The  noise  may  also  play  a  role  in  the  pitch-angle  scattering  and 
diffusion  of  energetic  electrons  observed  in  the  vicinity  of  the  ion 
cloud.  Because  of  the  similarity  to  the  solar  wind  interaction  with 
the  gaseous  envelope  of  a  comet,  it  is  expected  that  similar  types  of 
wave-particle  effects  may  occur  upstream  of  comets. 
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I.  INTRODUCTION 

On  September  11  and  again  on  September  20,  1984,  a  cloud  of 
lithium  gas  was  released  in  the  solar  wind  upstream  of  the  earth's 
bow  shock  as  part  of  the  AKPTE  (Active  Magnetospheric  JTracer 
Explorers)  program.  Numerous  effects  were  observed  in  the  immediate 
vicinity  of  the  releases,  including  (1)  the  creation  of  a  dense 

j 

rapidly  expanding  cloud  of  ionized  lithium,  (2)  the  formation  of  a 
diamagnetic  cavity,  and  (3)  the  occurrence  of  an  intense  shock-like 
burst  of  electrostatic  noise  in  the  region  upstream  of  the  cavity.  A 
summary  of  these  effects  is  given  in  an  accompanying  paper  by 
Haerendel  et  al.  [1985].  The  purpose  of  this  paper  is  to  present  a 
detailed  analysis  of  the  shock-like  electrostatic  noise  and  to 
investigate  the  origin  of  this  noise.  Comparisons  are  also  made  with 
a  similar  type  of  noise  observed  in  the  earth's  bow  shock. 

The  primary  objective  of  the  AMPTE  program  is  to  use  artifi¬ 
cially  injected  barium  and  lithium  ions  as  tracers  to  study  the 
transport  and  energization  of  charged  particles  in  the  earth's  mag¬ 
netosphere,  and  a  secondary  objective  is  to  study  the  physics  of 
the  injection  process  [Krimigis  et  al. ,  1982].  Three  spacecraft  are 
involved  in  the  project:  the  Ion  Release  Module  (IRM)  which  carries 
the  barium  and  lithium  canisters,  the  United  Kingdom  Subsatellite 
(UKS)  which  orbits  near  the  IRM  to  study  the  expansion  and  evolution 
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of  the  ion  cloud,  and  the  Charge  Composition  Explorer  (CCE)  which 
searches  for  the  tracer  ions  in  the  inner  regions  of  the  magneto¬ 
sphere. 

The  data  presented  in  this  paper  are  from  the  plasma  wave 
experiment  on  the  IRM.  The  IRM  is  in  a  low-inclination  highly 
eccentric  orbit  with  an  apogee  geocentric  radial  distance  of  19.0  Rg. 
During  the  September  1984  releases,  the  apogee  was  on  the  dayside  of 
the  earth  near  local  noon.  The  ion  cloud  is  produced  by  releasing 
two  canisters  containing  lithium  or  barium  in  diametrically  opposite 
directions  from  the  IRM.  Ten  minutes  after  the  release  the  two 
canisters  are  exploded  simultaneously  at  a  distance  of  about  1  km 
from  the  spacecraft.  The  atoms  vaporized  by  the  explosion  then  form 
an  expanding  cloud  of  neutral  gas  which  sweeps  over  the  spacecraft  a 
few  seconds  after  the  explosion.  Ultraviolet  radiation  from  the  sun 
subsequently  ionizes  the  neutral  gas,  thereby  producing  a  cloud  of 
ionized  gas.  For  a  discussion  of  the  physics  of  the  ion  cloud  forma¬ 
tion  and  expansion,  see  Haerendel  [1983,  1985]. 

A  variety  of  particle  and  field  instruments  on  the  IRM  are  used 
to  analyze  the  plasma  effects  associated  with  the  ionized  gas  cloud. 
For  a  list  of  these  instruments  and  their  characteristics,  see 
Krimigis  et  al.  [1982].  The  measurements  of  particular  interest  for 
this  paper  are  from  the  plasma  wave  instrument.  The  plasma  wave 
instrument  uses  a  47m  tip-to-tip  electric  dipole  antenna  for  electric 
field  measurements  and  two  search  coil  antennas  for  magnetic  field 


measurements.  The  signals  from  the  electric  and  magnetic  antennas 
are  processed  by  a  variety  of  receivers  and  spectrum  analyzers 
provided  by  three  groups:  the  University  of  Iowa,  the 
Max-Planck-Institut ,  and  the  Aerospace  Corporation.  For  a  complete 
description  of  this  instrumentation,  see  Hausler  et  al.  [1985a]. 


II.  CHARACTERISTICS  OF  THE  INTENSE  ELECTROSTATIC  NOISE 


The  overall  features  of  the  plasma  wave  observations  are  very 
similar  for  both  the  September  11th  and  September  20th  lithium 
releases.  Since  a  survey  of  the  plasma  wave  observations  is  given  in 
a  companion  paper  by  Hausler  et  al.  [1985b],  we  will  not  repeat  the 
detailed  description  of  both  events.  For  purposes  of  analyzing  the 
origin  of  the  intense  electrostatic  noise  it  is  sufficient  to  use  the 
September  20th  event  as  a  model.  The  plasma  interactions  during  this 
event  tend  to  be  the  least  complicated  and  most  readily  amenable  to  a 
theoretical  analysis. 

A  summary  of  the  plasma  density,  magnetic  field  and  plasma  wave 
electric  field  data  obtained  during  the  September  20th  event  is  shown 
in  Figure  1.  The  top  panel  shows  the  plasma  density  obtained  from  an 
analysis  of  electron  plasma  oscillations  observed  by  the  plasma  wave 
instrument  [Anderson  et  al. ,  1985],  the  second  panel  shows  the  mag¬ 
netic  field  from  the  magnetometer  [Luhr  et  al. ,  1985],  the  third 
panel  shows  the  electric  field  spectrum  from  the  wideband  receiver, 
and  the  bottom  panel  shows  the  electric  field  intensities  from  the 
16-channel  plasma  wave  spectrum  analyzer.  The  canisters  of  lithium 
were  exploded  at  0956:02  UT,  near  the  left  edge  of  the  plot.  About 
0.2  seconds  after  the  injection  the  plasma  density  jumps  up  to  about 
lO'5  cm--1,  the  magnetic  field  drops  to  near  zero,  and  an  impulsive 


disturbance  occurs  in  the  low  frequency  electric  field  channels. 

These  effects  mark  the  arrival  of  the  lithium  ion  cloud.  The 
depressed  magnetic  field  indicates  that  the  spacecraft  is  in  the 
diamagnetic  cavity  formed  by  the  dense  ion  cloud.  The  geometry  of 
this  cavity  is  illustrated  in  Figure  2.  As  the  lithium  ion  cloud 
expands  the  plasma  density  decreases.  The  cloud  also  starts  to  move 
downstream,  away  from  the  sun,  because  of  stresses  imposed  by  the 
solar  wind.  About  7  seconds  after  the  explosion  the  magnetic  field 
begins  to  increase.  At  this  time  the  cloud  has  moved  downstream  to 
the  point  where  the  spacecraft  is  passing  through  the  sunward 
boundary  of  the  diamagnetic  cavity.  The  spacecraft  then  enters  a 
region  of  enhanced  plasma  density  and  strong  magnetic  field.  The 
peak  strength  of  the  magnetic  field  is  about  six  times  stronger  than 
the  ambient  solar  wind  magnetic  field.  This  region  of  compressed 
plasma  and  strong  magnetic  field  is  caused  by  the  "pileup"  of  plasma 
around  the  nose  of  the  cloud  due  to  the  interaction  with  the  solar 
wind.  The  magnetic  field  is  believed  to  be  draped  around  the  diamag¬ 
netic  cavity  as  shown  in  Figure  2. 

As  the  ion  cloud  moves  farther  downstream,  a  region  of  very 
intense  electric  field  noise  is  encountered  about  17  seconds  after 
the  injection.  This  region  of  intense  electric  field  noise  is 
located  on  the  sunward  side  of  the  compression  region,  in  a  region 
where  the  magnetic  field  is  returning  to  the  pre-event  value.  The 
region  of  intense  electric  field  noise  also  appears  to  coincide  with 
a  brief  increase  in  the  plasma  density.  Because  considerable  uncer¬ 
tainty  exists  about  the  precise  shape  of  the  plasma  density  profile 


in  this  region  (see  Anderson  et  al.  [1985]),  this  portion  of  the 
plasma  density  profile  is  shown  by  dashed  lines  in  Figure  2.  The 
intense  electric  field  noise  continues  for  about  10  seconds,  slowly 
decreasing  in  intensity  as  the  magnetic  field  returns  to  the  upstream 
value. 

At  no  time  during  the  event  was  any  noise  comparable  to  that 
shown  in  Figure  1  detected  with  the  search  coil  magnetometers.  The 
absence  of  a  comparable  magnetic  component  indicates  that  the  noise 
is  electrostatic,  which  is  not  surprising  since  no  electromagnetic 
mode  exists  in  the  frequency  range  where  the  noise  is  observed 
(between  the  electron  cyclotron  frequency  and  the  electron  plasma 
frequency).  The  spectrum  of  the  electrostatic  noise  is  shown  in 
Figure  3.  Two  types  of  spectrums  are  shown,  labelled  peak  and 
average.  The  peak  is  the  maximum  electric  field  spectral  density 
observed  over  a  1.0  second  interval  and  the  average  is  the  electric 
field  spectral  density  averaged  over  a  50  msec  interval.  The  time 
for  the  spectrums  shown  in  Figure  3  is  0956:19.53  UT,  which  is  near 
the  point  of  maximum  intensity.  The  peak  and  average  broadband  field 
strengths,  integrated  from  31  Hz  to  178  kHz,  are  40. A  and  6.6  mVm-1, 
respectively.  These  intensities  are  among  the  most  intense  ever 
recorded  by  a  space  plasma  wave  experiment. 

The  noise  in  Figure  3  is  characterized  by  a  very  broad  peak 
extending  from  a  few  tens  of  Hz  up  to  a  few  hundred  kHz.  Above  a  few 
kHz  the  intensity  drops  very  rapidly,  to  near  the  instrument  noise 
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level  at.  frequencies  above  about  10  kHz.  The  peak  spectrum  shows  a 
broad  maximum  at  a  frequency  of  about  178  Hz,  and  the  average  spec¬ 
trum  shows  a  discontinuity  in  the  slope  at  about  1.78  kHz.  These  ’ 
detailed  tcatutes  are  probably  not  significant,  since  they  vary  some¬ 
what  from  event  to  event.  For  example,  the  corresponding  average 
spectrum  for  the  September  11th  event,  shown  in  Figure  4,  does  not 
show  a  discontinuity  in  the  slope.  Also,  the  average  spectrum  for 
the  September  11th  event  shows  a  peak  at  about  178  Hz,  whereas  the 
average  spectrum  for  the  September  20th  event  decreases  monotonically 
over  the  entire  frequency  range  measured. 

High  resolution  wideband  spectrums,  such  as  in  the  third  panel 
of  Figure  1,  show  that  the  noise  tends  to  be  somewhat  impulsive, 
consisting  of  many  short  bursts  lasting  only  a  fraction  of  a  second. 
This  impulsive  character  accounts  for  the  large  difference  between 
the  peak  and  average  spectrums  in  Figure  3.  The  impulsive  character 
of  the  noise  also  varies  somewhat  from  event  to  event.  The  September 
11th  event  in  Figure  4  has  a  smaller  difference  between  the  peak  and 
average  intensities,  which  shows  that  the  noise  was  less  impulsive 
during  this  event.  No  evidence  for  spin  modulation  is  evident  in  the 
electric  field  intensities,  which  indicates  that  the  electric  field, 
hence  wave  vector  direction,  tends  to  be  distributed  over  a  broad 
range  of  angles. 


V.  CONCLUSION 


We  have  presented  a  detailed  analysis  of  the  shock-like  electro¬ 
static  noise  observed  upstream  of  the  AMPTE  artificial  ion  clouds  and 
have  described  an  ion  beam-plasma  instability  that  can  account  for 
the  main  features  of  this  noise.  Comparisons  with  other  spacecraft 
measures  nts  show  that  this  noise  is  very  similar  to  the  electro¬ 
static  noise  observed  in  the  earth's  bow  shock.  This  similarity  sug¬ 
gests  that  a  shock  may  be  present  in  the  solar  wind  flow  upstream 
of  the  ion  cloud.  The  available  data  do  not  confirm  the  existence  of 
a  shock,  but  also  do  not  rule  out  the  possibility.  If  the  noise  is 
associated  with  a  shock,  then  the  shock  cannot  be  a  MHD  shock  since 
the  ion  cyclotron  radii  are  large  compared  to  the  radius  of  curvature 
of  the  discontinuity.  The  most  likely  possibility  is  that  the  noise 
is  associated  with  an  electrostatic  shock.  This  possibility  is 
supported  by  the  fact  that  laboratory  experiments  [Ikezi  et  al. , 

1973]  demonstrate  that  turbulent  electrostatic  shocks  can  be  produced 
by  injecting  an  ion  beam  into  a  plasma  at  rest,  a  physical  situation 
qualitatively  similar  to  injecting  an  ion  cloud  into  the  supersoni¬ 
cally  flowing  solar  wind. 

An  investigation  of  the  instabilities  caused  by  a  solar  wind 
proton  beam  streaming  through  a  nearly  stationary  lithium  ion  cloud 
shows  that  the  basic  features  of  the  electrostatic  noise  can  be 
accounted  for  oy  an  electrostatic  ion  beam-plasma  instability.  Using 
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shift  the  electron  drift  velocity  relative  to  that  given  by  Equation 
2,  which  assumes  zero  current.  Using  representative  values  for  the 
magnetic  field  variation,  AB  »  33  nT,  and  scale  size,  AL  =  200  km,  of 
the  magnetic  field  structure,  the  current  density  is  estimated  to  be 
0.14  pAm-^.  For  an  electron  density  of  8  cm- 3 ,  this  current  produces 
a  shift  in  the  electron  drift  velocity  of  about  25%.  This  shift  is 
an  upper  limit,  since  in  the  more  dense  region  of  the  cloud  the  shift 
required  to  produce  the  required  current  is  smaller.  Figure  13  shows 
the  effect  of  a  ±252  shift  in  the  electron  drift  velocity  for  the 
case  h'ec  =  Nli»  Somewhat  similar  results  are  obtained  for  the  case 
Nec  =  0.  As  can  be  seen,  a  positive  shift  in  the  drift  velocity 
(toward  Vgw)  increases  the  size  of  the  unstable  region,  and  a  nega¬ 
tive  shift  decreases  the  size  of  the  unstable  region.  Generally, 
the  effect  of  the  current  is  quite  small.  Therefore,  the  basic 
instability  is  an  ion  beam-driven  instability  and  not  a  current- 
driven  instability.  The  currents  that  exist  in  the  upstream  boundary 
region  are  simply  too  small  to  play  a  major  role  in  the  stability 


considerations. 


22 


Although  the  cold  electron  temperature  does  not  have  a  major 
effect  on  the  instability  boundaries,  it  does  affect  the  growth  rate 
somewhat.  The  maximum  growth  rate  is  shown  in  Figure  12  as  a  func¬ 
tion  of  the  density  ratio  N^^/Np.  The  solid  line  is  for  the  case 
Nec  =  0,  and  the  dashed  line  is  for  the  case  Nec  =  N^.  As  can  be 
seen  the  growth  rate  increases  when  the  cold  electrons  are  present. 
This  dependence  apparently  occurs  because  the  Landau  damping 
decreases  as  the  electron  temperature  decreases,  which  causes  higher 
growth  rates. 

It  can  be  easily  shown  that  the  growth  rates  in  Figure  12,  typi¬ 
cally  25  to  75  sec-1,  are  more  than  adequate  to  produce  the  observed 
intensities.  From  Figure  9  the  group  velocity,  3u>/9k  is  estimated  to 
range  from  about  5.3  to  50  km/sec.  Taking  a  representative  group 
velocity  of  50  km/sec  and  a  path  length  of  30  km  a  wave  packet  would 
propagate  through  the  unstable  region  in  0.6  sec.  For  a  representa¬ 
tive  growth  rate  of  25  sec-1,  the  wave  amplitude  would  grow  by  a 
factor  of  e15.  Waves  of  lower  group  velocity  grow  to  even  larger 
amplitudes.  These  simple  estimates  show  that  in  a  very  short  time 
the  waves  should  grow  to  very  large  amplitudes.  Although  the 
detailed  shapes  of  the  spectrums  in  Figures  3  and  4  are  controlled  by 
nonlinear  effects,  the  basic  frequency  range  where  the  noise  occurs 
appears  to  be  accounted  for  by  this  linear  stability  analysis. 

We  now  return  to  the  zero  current  assumption.  It  is  evident 
from  the  magnetic  field  variations  that  finite  currents  exist  in  the 
region  where  the  noise  occurs.  The  main  effect  of  a  current  is  to 
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A  representative  growth  rate  computation  for  the  instability 
between  0  and  is  shown  in  Figure  9.  The  bottom  panel  shows  the 
growth  rate,  y ,  plotted  as  a  function  of  frequency,  and  the  top  panel 
shows  the  wave  number,  kXjj,  normalized  by  the  Debye  length,  plotted 
as  a  function  of  frequency.  The  lithium  and  proton  density  are 
assumed  to  be  equal,  which  represents  the  most  unstable  condition. 

The  solid  curve  gives  the  case  where  Nec  =  0,  and  the  dashed  curve 
gives  the  case  where  Nec  =  N^.  Both  cases  give  large  growth  rates, 
25  sec-1  and  72  sec-1,  over  a  broad  range  of  frequencies,  up  to  about 
250  Hz.  This  range  of  frequencies  is  in  good  agreement  with  the 
spectrums  in  Figures  3  and  4,  which  show  that  the  most  intense  noise 
occurs  at  frequencies  from  a  few  tens  of  Hz  up  to  a  few  hundred  Hz. 

Because  the  relative  densities  of  the  lithium  ions  and  protons 
vary  over  a  large  range,  it  is  useful  to  investigate  the  instability 
as  a  function  of  the  lithium  to  proton  density  ratio,  Nj^/Np. 

Figures  10  and  11  show  the  frequency  of  marginal  stability  (y  -  0)  as 
a  function  of  the  density  ratio,  N^^/Np.  The  dashed  line  labelled 
Ymax  gives  the  frequency  of  maximum  growth  rate.  Figure  10  is  for 
the  case  Nec  =  0,  and  Figure  11  is  for  the  case  Nec  =  N^.  In  both 
cases  the  instability  occurs  over  a  broad  region,  which  shows  that 
the  basic  mechanism  is  relatively  insensitive  to  the  assumed  param¬ 
eters.  Figures  10  and  11  also  confirm  the  previously  mentioned  fact 
that  the  instability  only  occurs  over  a  limit  range  of  lithium  to 
proton  density  ratios.  If  the  lithium  to  proton  density  ratio  is 
less  than  about  2.5  *  10-^  the  instability  ceases,  and  if  the  ratio 
is  above  about  50  the  the  instability  again  ceases. 


dispersion  function  [Fried  and  Conte,  1961].  The  parameter  zs  is  a 
dimensionless  quantity  given  by 


z 

s 


(4) 


where  ms  is  the  mass  of  the  sth  component  and  k  is  Boltzmann's  con¬ 
stant.  The  dispersion  relation  has  been  solved  numerically  for  a)  and 
Y  using  the  Muller  method  [Muller,  1956].  Since  it  is  not  known 
whether  a  cold  electron  component  is  in  fact  present,  the  dispersion 
relation  was  analyzed  for  two  extreme  cases.  First,  it  was  assumed 
that  Nec  =  Nn*  This  case  represents  the  situation  where  cold  elec¬ 
trons  associated  with  the  lithium  photoionization  are  still  bound  to 
the  ion  cloud.  Second,  it  was  assumed  that  Nec  =  0.  This  amounts  to 
assuming  that  all  of  the  photoelectrons  have  escaped  from  the  cloud 
and  that  no  cold  electrons  are  present. 

Of  the  two  instabilities  discussed  earlier,  the  growth  rate 
computations  show  that  the  instability  with  phase  velocities  between 
0  and  V0  is  the  most  important.  The  reason  is  that  the  lithium  ion 
temperature,  TLj  ~  2  x  103°K,  is  much  lower  than  the  proton  tempera¬ 
ture,  Tp  =  1  x  105oK,  which  makes  the  double-hump  associated  with  the 
lithium  ions  and  electrons  much  more  pronounced  than  the  double  hump 
associated  with  the  electrons  and  protons.  The  instability  with 
phase  velocities  between  and  Vsw  only  occurs  for  proton  tempera¬ 
tures  less  than  about  3  x  10‘+oK.  Since  the  measured  proton  tempera¬ 
ture  is  1  x  105oK,  this  instability  never  occurs. 
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ambipolar  electric  field  eventually  disappears  and  the  photoelectrons 
escape.  The  photoelectron  density  is  then  zero,  Nec  =  0.  Macro¬ 
scopic  charge  neutrality  then  implies  that  Ne^  =  Np  +  NLi»  Finally, 
we  must  consider  the  electron  drift  velocity.  As  a  first  approxima- 
tion  the  E  *  B  drift  velocity  can  be  estimated  by  assuming  that  the 
net  current  in  the  plasma  is  zero.  Later  we  will  consider  the  effect 
of  a  finite  current.  Using  the  condition  for  macroscopic  charge 
neutrality,  the  zero  current  condition  gives 


(2) 


During  the  early  dense  phase  of  the  ion  cloud  expansion,  where  NL£  » 
Np,  the  electron  drift  velocity  is  near  zero,  and  during  the  later 
tenuous  phase,  where  N^i  <<  Np,  the  electron  drift  velocity  is  near 
the  solar  wind  velocity. 

Using  the  above  parameters,  the  frequency  w  and  growth  rate  y  of 
electrostatic  waves  can  be  calculated  from  the  electrostatic  disper¬ 
sion  relation,  which  for  a  Maxwellian  distribution  is 


D(  id,  y)  =  1  +  E  tt4 — rz  (1  +  zsZ(zs)]  =  0  ,  (3) 

s  ' K aDs ' 

where  the  summation  is  over  all  components  (s),  k  is  the  wave  number, 
aDs  is  the  Debye  length  of  the  sth  component,  and  Z  is  the  plasma 


proton  densities,  N^j[  and  Np,  which  are  shown  in  Figure  1,  (4)  the 
solar  wind  proton  temperature,  Tp,  which  is  about  105°K  [Coates  et 
al. ,  1985],  and  (5)  the  hot  electron  temperature,  Tgh,  which  at 
energies  above  about  50  eV  is  about  5  x  105°K  [Paschmann  et  al. , 
1985].  Parameters  that  are  unknown  because  they  cannot  be  measured 
by  the  plasma  instrument  are  (1)  the  lithium  ion  temperature,  TLi, 

(2)  the  cold  (photoelectron)  temperature,  Tec,  and  density  Nec,  and 

(3)  the  drift  velocity  of  the  electrons,  V^. 

To  proceed  further  suitable  choices  must  be  made  for  the  unknown 
parameters.  It  is  believed  that  the  chemical  reaction  responsible 
for  vaporizing  the  lithium  yields  a  temperature  of  about  2  x  103°K 
[Haerendel  et  al. ,  1985].  Although  the  lithium  ions  could  be  colder 
due  to  time-of-f light  effects,  we  will  adopt  the  value,  =  2  x 

103°K,  for  the  lithium  ion  temperature.  Since  nothing  is  known  about 
the  photoelectrons  associated  with  the  lithium  ionization,  the  best 
that  can  be  done  is  to  consider  two  extreme  situations.  During  the 
early  phase  of  the  ion  cloud  expansion  it  is  likely  that  the  photo¬ 
electrons  are  bound  to  the  ion  cloud  by  the  ambipolar  electric  field. 
In  this  case  the  cold  electron  density  is  equal  to  the  lithium  ion 
density,  Nec  =  N^.  Macroscopic  charge  neutrality  then  implies  that 
Ne^  =  Np.  Since  the  photoelectrons  are  ejected  from  the  lithium  ions 
with  an  energy  distribution  comparable  to  the  characteristic  energy 
of  the  solar  ultraviolet  spectrum  (~  1  eV)  it  is  likely  that  the  cold 
electron  temperature  is  about  104oK.  We  will  adopt  the  value  Tec  = 
104oK.  During  the  later  phase  of  the  ion  cloud  expansion  the 
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small,  V(j  is  expected  to  lie  between  the  drift  velocities  of  the 
lithium  ions  and  the  solar  wind  protons,  as  illustrated  in  Figure  8. 
The  reduced  distribution  function  then  has  two  double  humps,  one 
associated  with  the  lithium  ions  and  electrons,  and  the  other  asso¬ 
ciated  with  the  electrons  and  solar  wind  protons.  Two  instabilities 
are  then  expected,  one  with  phase  velocities  0  <  w/k  <  ,  and  the 

other  with  phase  velocities  Vd  <  co/k  <  Vsw.  These  instabilities  only 
occur  if  the  lithium  ion  and  solar  wind  proton  densities  are  somewhat 
comparable.  Otherwise,  the  double  hump  associated  with  the  minority 
ion  disappears  because  the  "beam"  density  is  too  low,  and  the  double¬ 
hump  associated  with  the  majority  ion  disappears  because  the  elec¬ 
trons  drift  at  essentially  the  same  velocity  as  the  majority  ion. 
These  basic  facts  provide  a  simple  explanation  for  why  the  noise  only 
occurs  near  the  outer  boundary  of  the  ion  cloud.  It  is  only  in  this 
region  that  the  lithium  ion  and  solar  wind  proton  densities  are 
comparable  (see  Figure  1). 

To  perform  a  detailed  stability  analysis  we  assume  that  the 
plasma  distribution  function  consists  of  four  Maxwellian  components 
representing  the  lithium  ions,  the  solar  wind  protons,  the  "hot" 
solar  wind  electrons,  and  the  "cold"  lithium  photoelectrons.  Each 
component  is  characterized  by  a  density,  a  temperature,  and  a  drift 
velocity.  For  the  September  20th  event  the  parameters  that  are 
reasonably  well  known  are  (1)  the  lithium  drift  velocity,  which  is 
essentially  zero,  (2)  the  solar  wind  velocity,  Vgw,  which  is  about 
460  km/sec  [Paschmann  et  al. ,  1985],  (3)  the  lithium  and  solar  wind 
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interaction,  first  analyzed  in  detail  by  Fried  and  Wong  [1966]. 

For  a  further  discussion  of  the  ion  beam-plasma  instability,  see  Imre 
and  Ozizmir  [1974]. 

It  can  be  easily  shown  that  the  growth  rate  of  electrostatic 
instabilities  in  an  unmagnetized  plasma  is  determined  by  a  reduced 
(one-dimensional)  distribution  function,  F(v),  which  is  the  sum  of 
the  electron  (e)  and  ion  (i)  distribution  functions, 

F(v)  =  Fe  +  I  —  Fi(v)  ,  (1) 

i  mi 

where  the  ion  terms  F^(v)  have  weighting  factors  given  by  the  elec¬ 
tron  to  ion  mass  ratio,  me/m^.  The  plasma  is  unstable  when  the 
reduced  distribution  function  is  sufficiently  double-humped  to 
satisfy  the  Penrose  criterion  [Krall  and  Trivelpiece,  1973].  Because 
the  electron  to  ion  mass  ratio  is  always  very  small,  the  ion  terms 
must  be  sharply  peaked  and  well  separated  for  instability  to  occur. 
Because  the  thermal  velocities  of  both  the  lithium  ions  and  the  solar 
wind  protons  are  much  smaller  than  the  solar  wind  velocity,  the 
distribution  functions  tend  to  satisfy  this  condition.  However,  it 
is  also  necessary  to  consider  the  role  of  the  electrons.  The 
electron  distribution  is  expected  to  consist  of  two  components:  a 
"cold"  component  associated  with  the  lithium  photoionization,  and  a 
"hot"  component  associated  with  the  solar  wind  electrons.  Because 
the  cyclotron  radius  of  the  electrons  is  small  compared  to  the 
diameter  of  the  cloud,  both  components  are  expected  to  drift  at  the 
E  x  b  drift  velocity,  Vj.  As  long  as  the  current  is  reasonably 


IV.  ORIGIN  OF  THE  ELECTROSTATIC  NOISE 


Since  the  solar  wind  protons  stream  almost  uninhibited  through 
the  ion  cloud,  it  is  evident  that  an  ion  beam-plasma  instability  may 
be  a  good  candidate  for  producing  the  electrostatic  noise.  This 
instability  would  then  be  similar  to  the  mechanism  responsible  for 
the  electrostatic  noise  in  the  earth's  bow  shock.  For  super-critical 
shocks,  such  as  in  Figure  5,  the  electrostatic  noise  is  closely 
correlated  with  the  onset  of  ions  (Nj)  reflected  from  the  bow  shock 
[Gurnett,  1985].  These  reflected  ions  produce  an  ion  beam-plasma 
interaction  that  is  very  similar  to  the  interaction  that  occurs  as 
the  solar  wind  protons  stream  through  a  stationary  ion  cloud. 

Before  we  consider  a  quantitative  analysis  of  the  ion  beam- 
plasma  instability  it  is  useful  to  discuss  the  qualitative  aspects  of 
this  instability.  The  ion  distribution  function  in  the  vicinity  of 
the  ion  cloud  is  illustrated  schematically  in  Figure  7.  The  lithium 
ions  have  a  drift  velocity  near  zero  since  they  are  born  essentially 
at  rest  in  the  spacecraft  frame  of  reference.  The  solar  wind  protons 
drift  at  a  velocity  corresponding  to  the  upstream  solar  wind  velo¬ 
city,  Vgw.  The  relative  density  of  the  lithium  ions  and  solar  wind 
protons  varies  with  position  in  the  ion  cloud.  Far  upstream  of  the 
cloud  the  solar  wind  protons  dominate,  and  in  the  dense  central 
region  of  the  cloud  the  lithium  ions  dominate.  The  double  ion 
distribution  in  Figure  7  corresponds  to  the  classical  ion  beam-plasma 
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cyclotron  radius  of  a  solar  wind  proton  is  56  km.  Lithium  ions 
injected  into  the  solar  wind  at  a  relative  velocity  of  400  km/sec 
have  an  even  larger  cyclotron  radius,  up  to  2,500  km.  Since  these 
ion  cyclotron  radii  are  larger  than  the  scale  size  of  the  ion  cloud, 
the  upstream  disturbance  cannot  be  a  MHD  shock. 

If  the  disturbance  is  not  a  MHD  shock  then  the  only  remaining 
possibility  is  that  it  is  an  electrostatic  shock.  Electrostatic 
shocks  are  characterized  by  a  jump  in  the  electrostatic  potential  on 
a  spatial  scale  small  compared  to  the  ion  cyclotron  radii.  Because 
the  characteristic  thickness  of  an  electrostatic  shock  can  be  as 
small  as  a  few  Debye  lengths,  which  is  only  about  one  hundred  meters 
in  the  solar  wind,  an  electrostatic  shock  could  in  principle  occur 
upstream  of  the  artificial  ion  cloud.  Electrostatic  shocks  have 
received  extensive  theoretical  study  [Montgomery  and  Joyce,  1969; 
Forslund  and  Shonk,  1970;  Tidman  and  Krall,  1971]  and  have  been 
produced  in  laboratory  plasmas  [Ikezi  et  al.,  1973;  Means  et  al. , 
1973].  Both  laminar  and  turbulent  electrostatic  shocks  have  been 
observed.  Turbulent  shocks  are  mainly  associated  with  high  Mach 
numbers,  in  which  case  the  turbulence  is  produced  by  ions  reflected 
from  the  electrostatic  potential  barrier  at  the  shock.  If  the  elec¬ 
trostatic  noise  observed  upstream  of  the  artificial  ion  cloud  is 
caused  by  an  electrostatic  shock,  it  is  almost  certainly  a  turbulent 
shock,  since  the  Mach  number  of  the  solar  wind  flow  (8  to  10)  is 


quite  high 


shows  that  the  solar  wind  protons  stream  through  the  ion  cloud. 
Second,  the  solar  wind  velocity  must  decrease  as  the  plasma  flows 
across  the  discontinuity.  This  condition  also  appears  to  be  satis¬ 
fied.  For  the  September  20th  event  the  UKS  plasma  instrument  shows 
that  the  solar  wind  velocity  decreases  by  about  50  km/sec  near  the 
upstream  edge  of  the  cloud  [Coates  et  al.  ,  1985],  The  September  11th 
event  shows  a  similar  decrease  in  the  solar  wind  velocity.  Third,  to 
qualify  as  a  shock  the  entropy  must  increase  as  the  plasma  flows 
across  the  discontinuity.  Wave-particle  interactions  caused  by  the 
intense  electric  field  noise  upstream  of  the  ion  cloud  certainly 
could  cause  an  increase  in  entropy.  However,  the  entropy  production 
cannot  be  verified  directly  because  the  plasma  instrument  does  not 
have  sufficient  space  and  time  resolution  to  accurately  determine  the 
temperature  and  density  variations  in  the  region  where  the  noise  is 
observed. 

The  above  considerations  support  the  view  that  the  noise  up¬ 
stream  of  the  ion  cloud  could  be  associated  with  a  shock,  but  do  not 
prove  that  it  is  a  shock.  If  it  is  a  shock,  then  it  cannot  be  a  con¬ 
ventional  magnetohydrodynamic  (MHD)  shock.  For  a  MHD  shock  the  ion 
cyclotron  radii  must  be  small  compared  to  the  radius  of  curvature  of 
the  discontinuity,  otherwise  magnetohydrodynamics  would  not  apply. 
From  relatively  simple  pressure  balance  considerations  Haerendel  et 
al.  [1985]  estimate  that  the  radius  of  the  diamagnetic  cavity  is 
about  30  km.  For  a  representative  solar  wind  ion  temperature  of 
105oK  and  a  magnetic  field  strength  of  8nT  it  can  be  shown  that  the 


quite  striking.  In  both  cases  the  noise  is  almost  completely  elec¬ 
trostatic  and  extends  over  a  broad  range  of  frequencies,  from  a  few 
Hz  to  several  kHz.  An  electric  field  spectrum  of  the  noise  observed 
in  the  earth's  bow  shock  is  shown  Figure  6.  This  spectrum  is  from 
the  time  of  maximum  intensity  at  0251:20  UT.  The  similarity  of  this 
spectrum  to  the  spectrums  in  Figures  3  and  4  is  clearly  evident.  In 
both  cases  the  noise  extends  with  relatively  high  intensities  from  a 
few  Hz  to  a  few  kHz,  and  then  decreases  rapidly  with  increasing  fre¬ 
quency.  As  in  the  case  of  the  noise  upstream  of  the  artificial  ion 
cloud,  the  detailed  shape  of  the  spectrum  in  the  bow  shock  varies 
somewhat  from  event  to  event.  In  some  cases  the  spectrum  has  a  peak 
at  a  frequency  of  a  few  hundred  Hz,  similar  to  the  example  in  Figure 
4,  whereas  in  other  cases  it  is  nearly  flat,  as  in  Figure  6.  Wide¬ 
band  spectrums  show  that  the  noise  in  the  bow  shock  is  very  impulsive 
and  usually  does  not  have  a  detectable  spin  modulation,  very  similar 
to  the  noise  observed  upstream  of  the  artificial  ion  cloud.  For  a 
survey  of  representative  spectrums  of  the  electric  field  noise  in  the 
earth's  bow  shock,  see  Rodriguez  and  Gurnett  [1975].  Although  the 
overall  features  are  similar,  the  noise  associated  with  the  artifi¬ 
cial  ion  cloud  is  more  intense  than  the  noise  in  the  bow  shock, 
typically  by  a  factor  of  ten  to  one  hundred. 

To  interpret  the  effects  upstream  of  the  ion  cloud  as  a  shock 
certain  basic  conditions  must  be  satisfied.  First,  the  plasma  must 
have  a  component  of  flow  normal  to  the  surface  of  the  discontinuity. 
This  condition  is  known  to  be  satisfied  since  the  plasma  instrument 
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III.  IS  THE  ELECTROSTATIC  NOISE  ASSOCIATED  WITH  A  SHOCK? 

The  location  of  the  intense  electrostatic  noise  upstream  of  the 
ion  cloud  raises  the  question  as  to  whether  this  noise  is  associated 
with  a  shock.  For  many  years  it  has  been  known  that  an  intense  burst 
of  electrostatic  noise  is  present  in  the  earth's  bow  shock  [Fredricks 
et  al.  ,  1968,  1970a, b;  Rodriguez  and  Gurnett,  1975].  It  is  now 
widely  believed  that  in  a  collisionless  plasma  this  noise  causes 
heating  and  dissipation  very  similar  to  collisions  in  an  ordinary 
gas-dynamic  shock.  The  close  similarity  of  the  electrostatic  noise 
to  the  noise  in  the  earth's  bow  shock  is  illustrated  in  Figure  5, 
which  shows  a  crossing  of  the  earth's  bow  shock  by  the  ISEE  1  space¬ 
craft  on  November  7,  1977.  The  top  panel  of  Figure  5  shows  the  elec¬ 
tric  field  intensities  and  the  bottom  three  panels,  from  Paschmann  et 
al.  [1982],  show  the  electron  density  Ne,  the  reflected  ion  density, 
Nj,  and  the  magnetic  field  strength,  B.  The  shock  can  be  identified 
by  the  abrupt  increase  in  the  electron  density  and  magnetic  field 
strength  from  about  2250:55  to  2251:20  UT.  This  region  of  rapidly 
increasing  plasma  density  and  magnetic  field  strength  is  called  the 
transition  region.  The  top  panel  of  Figure  5  shows  that  an  intense 
burst  of  electric  field  noise  is  present  in  the  transition  region. 

The  similarities  between  the  electric  field  noise  upstream  of 
the  artificial  ion  cloud  and  the  noise  in  the  earth's  bow  shock  are 


realistic  parameters  it  is  found  that  the  ion  beam-plasma  instability 
occurs  for  a  range  of  lithium-to-proton  density  ratios  between  about 
2.5  x  10-2  and  50.  These  limits  explain  why  the  noise  only  occurs 
near  the  outer  edge  of  the  ion  cloud,  in  the  region  where  the 
lithium-to-proton  density  ratio  is  near  unity.  The  maximum  growth 
rates  increase  if  cold  electrons  associated  with  the  lithium  photo¬ 
ionization  are  present  in  the  plasma.  Currents  do  not  play  an 
important  role  in  generating  the  noise. 

At  present  the  importance  of  the  electrostatic  noise  has  not 
been  established.  The  noise  is  very  intense  and  occurs  in  a  region 
of  elevated  electron  temperatures.  Although  the  main  electron 
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energization  is  thought  to  be  caused  by  adiabatic  acceleration  in  the 
magnetic  field  upstream  of  the  cloud  [Paschmann  et  al. ,  1985],  the 
noise  may  play  a  role  in  the  pitch-angle  scattering  and  diffusion  of 
these  electrons.  If  the  noise  is  associated  with  an  electrostatic 
shock,  then  the  turbulent  electric  fields  probably  play  an  important 
role  in  the  ion  heating  and  dissipation  in  the  shock.  Since  the 
phase  velocity  of  the  unstable  waves  is  near  the  lithium  thermal 
velocity,  the  turbulence  should  be  very  effective  at  heating  the 
lithium  ions. 

It  is  Likely  that  the  intense  electrostatic  noise  observed  dur¬ 
ing  the  AMPTE  ion  releases  is  closely  related  to  a  variety  of  other 
situations  where  intense  electric  field  noise  is  produced  by  a 
neutral  gas  interacting  with  a  rapidly  moving  plasma.  These  situa¬ 
tions  include  (1)  intense  electric  field  noise  associated  with  cesium 
ion  release  in  the  earth's  ionosphere  [Kintner  et  al.,  I980J, 


(2)  intense  electrostatic  noise  observed  in  the  vicinity  of  the  Space 
Shuttle  [Shawhan  and  Murphy,  1983],  apparently  caused  by  an  inter¬ 
action  of  gases  from  the  Shuttle  with  the  surrounding  ionosphere,  and 

(3)  electric  field  noise  observed  by  Voyager  1  in  the  vicinity  of 
Titan  [Gurnett  et  al. ,  1981].  Because  of  the  close  similarity  of  the 
lithium  releases  to  the  gaseous  envelope  around  a  comet,  it  will  be 
interesting  to  see  if  a  comparable  type  of  noise  is  observed  during 
the  forthcoming  fly-bys  of  the  comets  Giacobini-Zinner  and  Halley. 
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FIGURE  CAPTIONS 


Figure  1  A  summary  of  Che  plasma  density,  magnetic  field, 

and  electric  field  observations  for  the  September 
20,  1984,  solar  wind  lithium  release.  Note  the 
intense  burst  of  electrostatic  noise  from  about  19 
to  28  seconds  after  the  injection,  as  the  magnetic 
field  returns  to  the  upstream  value.  A  somewhat 
similar  burst  of  noise  was  observed  following  the 
September  11,  1984,  lithium  release. 

Figure  2  A  schematic  illustration  showing  the  inferred 

plasma  and  magnetic  field  configuration  associated 
with  the  September  solar  wind  ion  releases.  The 
intense  burst  of  electrostatic  noise  is  observed  on 
the  upstream  edge  of  the  region  of  draped  magnetic 
field  lines. 

Figure  3  The  peak  and  average  electric  field  spectrums  for 

the  September  20th  event.  The  peak  spectrum  shows 
that  the  maximum  electric  field  intensities  occur¬ 
red  in  the  frequency  range  around  a  few  hundred  Hz. 
The  intensities  decrease  rapidly  at  frequencies 


above  a  few  kHz 
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Figure  4 


Figure  5 


Figure  6 


Figure  7 
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The  peak  and  average  electric  field  spectruras  for 
the  September  11th  event.  Again,  the  peak  electric 
field  intensities  are  in  the  region  around  one 
hundred  Hz,  with  significant  noise  intensities 
extending  up  to  a  few  kHz. 

The  plasma  density,  magnetic  field,  and  electric 
field  measurements  for  a  crossing  of  the  earth's  bow 
shock  by  the  ISEE-1  spacecraft  on  November  7,  1977. 
The  electrostatic  noise  observed  on  this  and  other 
shock  crossings  is  very  similar  to  the  noise  observed 
upstream  of  the  AMPTE  artificial  ion  clouds. 

The  electric  field  spectrum  at  the  time  of  maximum 
intensity  for  the  shock  crossing  in  Figure  5.  Note 
the  very  close  similarity  to  the  artificial  ion  cloud 
spectruras  in  Figures  3  and  4. 

A  schematic  representation  of  the  ion  distribution 
function  in  the  vicinity  of  the  artificial  ion  cloud. 
The  lithium  ions  are  born  essentially  at  rest  in  the 
spacecraft  frame  of  reference.  The  solar  wind 
protons  stream  essentially  uninhibited  through  the 
ion  cloud  with  a  velocity  Vsw.  This  double-humped 


distribution  is  known  to  be  highly  unstable 


Figure  8 


A  model  for  Che  cold  electron  distribution  in 


which  the  cold  and  hot  electrons  both  drift  at  a 
velocity,  Vq,  that  is  determined  by  the  v  x  B 
electric  field.  This  model  applies  mainly  to  the 
motion  perpendicular  to  the  magnetic  field. 

Figure  9  Plots  of  the  growth  rate,  y,  and  normalized  wave 

number,  kXp,  as  a  function  of  frequency  for  the 
instability  associated  with  the  minimum  of  F(v) 
between  0  and  Vq.  The  solid  lines  are  for  Nec  =  0, 
and  the  dashed  lines  are  for  Nec  = 

Figure  10  The  frequencies  of  marginal  stability,  y  =  0,  and 

maximum  growth  rate,  Ymax,  as  a  function  of  the 
lithium  to  proton  density  ratio,  N^/Np.  This  case 
is  for  Nec  =  0.  Note  the  limited  range  of  N^^/Np 
over  which  instability  occurs.  Also,  note  that  the 
maximum  frequency  of  the  instability  compares  favor¬ 
ably  with  the  observed  frequency  range  of  the  noise 
in  Figures  3  and  4. 

Figure  LI  The  frequencies  of  marginal  stability,  y  =  0,  and 

maximum  growth  rate,  ymax,  as  a  function  of  the 
lithium  to  proton  density  ratio,  N^^/Np.  This  case 
is  for  Nec  = 
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Figure  12 


Figure  13 


The  maximum  growth  rate  of  the  ion  beam-plasma 
instability  as  a  function  of  the  lithium  to  proton 
density  ratio.  Note  that  the  highest  growth  rate 
occurs  when  Nj^  =  Np.  This  result  also  agrees  with 
the  observations  in  Figure  1,  which  show  that  the 
most  intense  noise  occurs  at  about  the  time  that 

NLi  “  V 

Marginal  instability  boundaries  comparable  to 
Figure  13  showing  the  effect  of  a  current.  The 
cross-hatched  region  shows  the  variation  in  the  y  = 
0  boundary  for  a  ±25%  variation  in  the  electron 
drift  velocity.  These  limits  represent  the 


maximum  variations  that  could  be  produced  by  the 
currents  in  the  upstream  boundary  region. 
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